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Dissociative  adsorption  of  molecular  oxygen  on  the  Al(lll)  surface  exhibits  mechanistic  complexity 
that  remains  surprisingly  poorly  understood  in  terms  of  the  underlying  physics.  Experiments  clearly 
indicate  substantial  energy  barriers  and  a  mysteriously  large  number  of  adsorbed  single  oxygen  atoms 
instead  of  pairs.  Conventional  first  principles  quantum  mechanics  (density  functional  theory)  predicts  no 
energy  barrier  at  all;  instead,  spin  selection  rules  have  been  invoked  to  explain  the  barrier.  In  this  Letter, 
we  show  that  correct  barriers  arise  naturally  when  embedded  correlated  electron  wave  functions  are  used 
to  capture  the  physics  of  the  interaction  of  02  with  the  metal  surface.  The  barrier  originates  from  an  abrupt 
charge  transfer  (from  metal  to  oxygen),  which  is  properly  treated  within  correlated  wave  function  theory 
but  not  within  conventional  density  functional  theory.  Our  potential  energy  surfaces  also  identify  oxygen 
atom  abstraction  as  the  dominant  reaction  pathway  at  low  incident  energies,  consistent  with  measure¬ 
ments,  and  show  that  charge  transfer  occurs  in  a  stepwise  fashion. 
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Understanding  the  interaction  of  gas  molecules  with 
metal  surfaces  is  critically  important  for  many  applications 
including  molecular  sensing,  catalysis,  and  corrosion.  The 
initial  dissociative  sticking  of  gas  molecules  such  as  mo¬ 
lecular  oxygen  (02)  enables  subsequent  reactions  through 
catalysis  [1,2].  Conversely,  surface  corrosion  engenders 
considerable  costs  [3].  Despite  their  importance,  the  details 
of  these  processes  are  poorly  understood  [4].  Even  for  the 
simple  yet  prominent  benchmark  case  of  02  approaching 
an  ideal  aluminum  (111)  surface  [Fig.  1(a)],  there  is  a 
troubling  discrepancy  between  theory  and  experiment: 
experiments  consistently  suggest  the  existence  of  an 
energy  barrier  [5-7]  preventing  low-energy  sticking. 
However,  studies  [8-10]  employing  conventional  Kohn- 
Sham  density  functional  theory  (KS  DFT)  show  no  such 
barrier.  There  is  an  ongoing  debate  whether  the  measured 
energy  barrier  is  due  to  spin  quenching  or  charge  transfer 
[4].  Furthermore,  experiments  consistently  find  widely 
spaced  oxygen  atoms  instead  of  pairs  (as  one  would  expect 
from  full  chemisorption  of  02),  especially  at  low  incident 
energies  [6,11].  Various  mechanisms  to  explain  these 
observations  have  been  proposed,  including  abstraction 
[6,12],  dissociation  [13],  and  hot-atom  motion  [11,14]. 

Recent  constrained  KS  DFT  calculations  related  the  ori¬ 
gin  of  the  energy  barrier  to  spin  selection  rules  [15,16],  as 
suggested  by  Ref.  [7].  By  contrast,  earlier  simulations 
employing  semiempirical  potentials  proposed  nonadiabatic 
charge  transfer  from  the  surface  to  02  [17].  Subsequent 
calculations  of  the  interaction  between  02  and  small  A1 
clusters  support  the  latter  model:  using  correlated  wave 
function  methods,  they  suggest  the  approximate  treatment 
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of  exchange  correlation  (XC)  in  KS  DFT  as  reason  for  its 
failure  [18-20].  Local  (or  semilocal)  density-based  XC 
functionals  lack  derivative  discontinuities  [21]  and  suffer 
from  self-interaction  error  due  to  their  lack  of  exact  ex¬ 
change  [22] .  These  two  XC  flaws  within  pure  DFT  produce 
overdelocalization  of  electrons.  We  therefore  expect  an 
unphysically  facile  charge  transfer  between  02  and  the 
A1  surface  with  conventional  KS  DFT  [18-20].  To  obtain 
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FIG.  1  (color  online),  (a)  Four-layer  Al(lll)  periodic  slab  and 
approaching  02  molecule,  (b)  Geometries  used  in  this  work  to 
investigate  different  points  of  the  Al(lll)  surface.  Both  perpen¬ 
dicular  and  parallel  incidence  of  02  were  investigated, 
(c-f)  The  four  different  high-symmetry  positions  of  the 
Al(lll)  surface,  and  the  respective  embedded  two-layer  clusters 
used  in  our  calculations. 
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a  reliable  picture  of  the  charge  transfer  process,  one  would 
need  to  apply  correlated  wave  function  methods  to  the 
entire  02/Al  system  [Fig.  1(a)],  which  is  too  computation¬ 
ally  demanding  to  be  feasible.  In  this  Letter,  we  solve  this 
problem  by  using  a  recently  developed  embedding  theory 
[23]  that  allows  for  a  completely  ab  initio  description  of 
the  interaction  between  02  and  its  nearby  A1  atoms  at  the 
different  possible  adsorption  sites  [Figs.  l(c)-l(f)],  while 
the  rest  of  the  A1  surface  is  treated  by  KS  DFT.  The  two 
subsystems  interact  via  a  unique,  formally  exact  embed¬ 
ding  potential  [23].  We  aim  for  a  detailed  understanding  of 
the  reaction,  including  possible  reaction  pathways  and 
their  associated  energy  barriers.  Our  results  provide  novel 
quantitative  insights  into  previous  measurements  [5,6]  and 
implicate  activated  charge  transfer  as  fundamental  to  the 
02/Al  adsorption  process.  The  highly  accurate  ab  initio 
potential  energy  surfaces  obtained  in  this  work  may  also 
serve  as  a  basis  for  future  dynamics  simulations  of  the 
reaction  of  oxygen  with  aluminum  surfaces. 

We  consider  an  02  molecule  approaching  an  Al(lll) 
surface  at  different  orientations  and  adsorption  sites 
[Fig.  1].  We  map  out  the  ground  state  energy  as  a  function 
of  the  02  bond  length  L0- o  and  the  distance  of  the  mole¬ 
cule’s  center-of-mass  to  the  A1  surface  LA1_o2  [see,  e.g., 
Fig.  2(a)].  For  a  pure  KS  DFT  description,  we  use  an  A1 
slab  [Fig.  1(a)]  described  by  a  periodic  5X5  supercell  (for 
details  see  Refs.  [24-26]).  To  elucidate  the  role  of  ex¬ 
change  and  correlation,  we  then  treat  the  same  problem 
with  highly  accurate  correlated  wave  function  (CW)  meth¬ 
ods  (second-order  multireference  many-body  perturbation 
theory  [26,27]).  We  cut  a  small  cluster  of  10  to  14  atoms 
out  of  the  A1  slab  to  model  the  different  possible  adsorption 


sites  [Figs.  l(c)-l(f)].  We  then  calculate  an  embedding 
potential  Vemb  representing  the  interaction  between  the 
cluster  and  the  remaining  Al(l  11)  surface  [23].  The  ground 
state  energy  of  this  embedded  cluster  in  the  presence  of  02 
is  then  determined  both  on  a  KS  DFT  (£^b)  as  well  as  a 
CW  (E^)  level.  Comparing  these  two  energies  yields  a 
correction  accounting  for  XC  effects  beyond  the  pure  KS 
DFT  slab  calculation  EDFT ,  so  that  the  final  embedded 
energy  £cmb  can  be  written  as  [23] 

F  =  FDFT  _i_  rCW  _  Z7 DFT  /i\ 

^emb  ^  “  ^emb  ^emb*  t 1 ) 

Several  software  packages  [24,27-32]  were  used  to  evalu¬ 
ate  Eq.  (1).  For  technical  details,  we  refer  the  interested 
reader  to  Ref.  [23]  and  the  supplemental  material  Ref.  [26]. 

As  discussed  in  the  introduction,  conventional  KS  DFT 
calculations  with  local  (or  semilocal)  density-based  XC 
functionals  cannot  correctly  describe  the  charge  transfer 
process  involved  in  chemisorption  of  02  on  Al(lll). 
Accordingly,  we  find  no  energy  barrier  in  our  KS  DFT 
slab  calculations  [Fig.  2(a)],  and  an  unphysically  smooth 
transfer  of  charge  from  the  metal  surface  to  the  02  mole¬ 
cule  (not  shown).  Our  results  change  qualitatively  upon 
considering  our  embedded  CW  energy  [Eq.  (1)]:  the  result¬ 
ing  embedding-based  potential  energy  surfaces  (EPESs) 
feature  pronounced  energy  barriers  for  all  surface  sites  we 
consider  [Figs.  2(b)  and  2(c)  and  Table  I].  The  maximum 
cluster-size-converged  barrier  height  (660  meV,  Table  I) 
appears  consistent  with  experiments  measuring  the  initial 
sticking  probability  of  02  on  Al(lll)  as  function  of  02 
kinetic  energy.  Almost  no  increase  of  sticking  probability 
is  observed  at  energies  above  600  meV  [5].  Fluctuations  in 
the  CW  results  due  to  inaccuracies  of  the  perturbative 
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FIG.  2  (color  online),  (a)  PES  calculated  using  KS  DFT  with  a  periodic  slab  model  for  02  approaching  perpendicular  to  the  fee 
hollow  site  on  the  Al(lll)  surface  [see  inset  and  Fig.  1].  (b)  Same  02  site  and  orientation  as  (a),  now  based  on  embedded  correlated 
wave  functions  [Eq.  (1),  6/6  cluster  geometry,  see  inset].  Barrier  height  is  360  meV  (Table  I).  Black  triangle  marks  an  intermediate 
local  minimum  and  dashed  line  an  elongated  minimum  within  a  reaction  pathway  towards  02  abstraction,  (c)  Same  method  as  (b)  for 
the  Al(lll)  bridge  site  with  02  parallel  to  the  surface  (8/4  cluster  geometry,  see  inset).  Barrier  height  is  530  meV  (see  Table  I  and 
Fig.  1  in  Ref.  [26]).  Increasing  cluster  size  to  8/6  increases  the  barrier  height  by  only  30  meV  (see  Table  I).  Minimum  of  the 
embedding -based  PES  at  LA1_0o  =  0.9 A  and  Lq-o  =  2.2  A.  Energy  spacing  between  contour  lines  is  300  meV. 
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TABLE  I.  Embedded  CW  predictions  of  the  barrier  height  and  surface  distance  at  the  top  of 
the  barrier  for  different  adsorption  sites  (Fig.  1),  cluster  sizes,  and  02  orientations.  Errors 
represent  fluctuations  in  the  final  EPES. 


Site 

No.  of  atoms 

1st  layer  2nd  layer 

Orientation 

^barrier  ±  $  (meV) 

(A) 

bridge 

8 

6 

ii 

560  ±  10 

2.4 

8 

6 

i 

450  ±  30 

2.7 

fee 

6 

6 

ii 

430  ±  20 

2.4 

6 

6 

i 

360  ±  10 

2.9 

hep 

6 

6 

ii 

410  ±  10 

2.5 

6 

6 

i 

410  ±  20 

2.9 

7 

3 

ii 

660  ±  20 

2.6 

top 

experiment  [5] 

7 

3 

i 

660  ±  20 
<  600 

2.8 

expansion  are  below  30  meV  (Table  I).  We  note  that  too- 
small  clusters  (with  less  than  ten  atoms)  overestimate  the 
energy  barrier  for  the  charge  transfer  [26],  due  to  the  smaller 
amount  of  polarizable  screening  (Vemb  is  not  self- 
consistently  updated  during  charge  transfer).  Calculations 
of  clusters  larger  than  14  atoms  quickly  become  unfeasible 
due  to  the  computational  demands  of  the  complete  active 
space  self-consistent  field  method  [26].  We  also  calculate 
energy  barriers  for  several  different  cluster  geometries  at 
the  same  adsorption  site  and  find  variations  ^  60  meV 
[26].  We  thus  estimate  a  final  uncertainty  of  ~  100  meV 
in  the  barrier  heights  caused  by  the  finite  cluster  size  and 
the  nonself-consistent  Vemb  [26]. 

To  investigate  the  origin  of  the  energy  barrier,  we  cal¬ 
culate  the  charge  and  spin  states  of  the  oxygen  molecule  as 
it  approaches  the  surface.  We  find  an  abrupt  charge  trans¬ 
fer,  accompanied  by  a  change  of  02  spin,  at  the  barrier 
[Fig.  3];  i.e.,  the  corresponding  diabatic  surfaces  cross. 
Far  from  the  A1  surface,  the  triplet  state  is  energetically 
favored.  At  LA1_o2  ~  2.5  A  [Table  I],  a  first  charge  transfer 
occurs  to  form  02~  and  a  doublet  spin  configuration  at 
the  02  becomes  preferred  [Fig.  3].  The  transfer  starts  at 
larger  LM-0l  for  larger  L0_0,  consistent  with  the  increased 
sticking  probability  of  vibrationally  excited  02  found  in 
experiment  [5]. 

To  assess  the  importance  of  including  an  embedding 
potential,  we  compare  (see  Table  1  in  Supplemental 
Material  Ref.  [26])  the  EPES  with  data  obtained  from  an 
isolated,  nonembedded  A1  cluster  treated  at  the  CW  level. 
While  abrupt  charge  transfer  at  the  crossing  of  the  barrier  is 
still  predicted,  the  relative  barrier  heights  are  incorrect:  the 
barrier  for  abstraction  is  much  higher  than  for  dissociative 
adsorption,  which  is  not  consistent  with  (and  is  in  fact  the 
reverse  of)  experiment.  The  embedding  potential  provides 
the  correct  boundary  condition  at  the  edges  of  the  cluster 
needed  to  obtain  the  correct  electronic  structure  inside  the 
A1  cluster.  The  correct  boundary  condition  is  key  to  obtain¬ 
ing  correct  trends  and  quantitative  results. 


After  crossing  the  barrier,  one  or  both  of  the  oxygen  atoms 
may  adsorb.  To  investigate  possible  reaction  pathways,  we 
calculate  the  energy  E0  [using  Eq.  (1)]  of  a  single  oxygen 
atom  at  various  surface  positions.  We  find  two  minima  at  the 
fee  (hep)  positions,  separated  by  1.6  A  and  with  energies  of 
8.3  eV  (8.1  eV)  lower  than  the  energy  of  an  oxygen  atom  at 
5  A  above  the  surface.  For  02  parallel  incidence  upon  the 
bridge  site,  both  oxygen  atoms  will  be  aligned  with  these 
minima.  Indeed,  the  corresponding  EPES  for  02  dissociation 
features  an  associated  minimum,  albeit  with  a  smaller 
^emb  =  —  6.8  eV  at  a  larger  L0-o  =  2.2  A  [Fig.  2(c)],  due 
to  Coulomb  repulsion  between  the  two  charged  oxygen 
atoms  [26]. 

The  above  energy  considerations  allow  for  a  second 
possible  reaction  pathway:  the  energy  gain  of  8.3  eV 
upon  adsorption  of  one  O  atom  is  large  enough  to  permit 
dissociative  adsorption  of  one  atom  while  the  second  atom 
is  emitted  into  the  gas  phase  as  a  free  radical,  a  process 
commonly  referred  to  as  abstraction.  Let  us  consider  the 
fee  site  and  perpendicular  incidence,  as  the  closer  oxygen 
atom  will  be  in  a  favorable  adsorption  position.  The  cor¬ 
responding  EPES  indeed  suggests  abstraction  for 
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FIG.  3  (color  online),  (a)  Mulliken  charge  populations  and 
(b)  Mulliken  spin  populations  for  02  approaching  the  bridge 
site  in  parallel  incidence.  The  contour  lines  of  the  potential 
energy  surface  [Fig.  2(c)]  are  superposed  to  guide  the  eye. 
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perpendicular  incidence:  we  observe  a  minimum  extending 
up  to  L0- o  =  2.5  A  [see  dashed  line  in  Fig.  2(b)]  that  is  not 
present  for  KS  DFT  [Fig.  2(a)].  The  relative  barrier  heights 
of  perpendicular  and  parallel  incidence  for  the  embedded 
CW  calculations  [Table  I]  suggest  abstraction  is  likely  to 
be  a  dominant  pathway  at  low  energies,  consistent  with  the 
observation  of  a  large  fraction  of  single  oxygen  atoms  on 
the  metal  surface  in  STM  measurements  at  low  incident 
energies  [6].  As  in  the  parallel  case,  we  find  abrupt  transfer 
of  charge  and  change  of  spin  [Fig.  4].  However,  now  there  is 
a  strong  asymmetry  between  charges  on  the  two  O  atoms: 
only  small  amounts  of  charge  are  transferred  to  the  farther 
oxygen  [Fig.  4(b)]  that  returns  to  a  triplet  configuration  at 
L0-o  >1.6  A  [Fig.  4(d)],  in  line  with  the 3 P  ground  state  of 
an  isolated  oxygen  atom,  leaving  the  adsorbed  oxygen  and 
the  cluster  in  a  spin-free  singlet  state.  Conversely,  the  closer 
oxygen  takes  nearly  the  same  charge  as  the  entire  02  mole¬ 
cule  in  the  parallel  case  [compare  Figs.  3(a)  and  4(a)].  Thus, 
the  EPES  suggests  that  abstraction  happens  in  a  three- 
step  process:  (i)  crossing  the  initial  barrier,  0.5  e  charge 
transfers  to  the  closer  atom;  (ii)  another  0.5  e  transfers 
to  form  C>2  in  a  shallow  local  minimum  [17]  of  the  potential 
energy  surface  [see  triangle  in  Fig.  4];  and  finally, 
(iii)  dissociation  of  a  neutral  O  in  its  3  P  ground  state,  while 
the  remaining  O  bonds  with  the  A1  surface  as  a  negatively 
charged  adsorbate. 

Finally,  we  compare  our  results  to  previous  theoretical 
models,  in  particular  nonadiabatic  [17]  or  diabatic  [16] 
approaches,  which  also  find  an  energy  barrier  consistent 
with  experiment.  However,  unlike  the  earlier  reports,  our 
description  contains  no  adjustable  model  parameters  or 
specially  designed  XC  functionals;  the  barriers  arise 


FIG.  4  (color  online).  Mulliken  charge  [(a-b)]  and  spin  [(c-d)] 
populations  for  close  and  far  O  atoms  in  02  in  perpendicular 
incidence  above  the  fee  site  [Fig.  2(b),  contour  lines  of  the  EPES 
with  a  line  spacing  of  400  meV  are  plotted  in  black  to  guide  the 
eye].  Black  triangles  mark  an  intermediate  local  minimum  along 
a  reaction  pathway  towards  02  abstraction.  Colors  refer  to 
legend  given  in  Fig.  3. 


naturally  when  electron  exchange  and  correlation  are  prop¬ 
erly  accounted  for.  Secondly,  our  relative  barrier  heights 
are  consistent  both  with  the  experimental  observation  of 
singly  absorbed  oxygen  atoms  at  low  incident  kinetic 
energies  and  the  observed  maximum  energy  threshold  for 
sticking.  By  contrast,  spin-restricted  calculations  overesti¬ 
mate  the  barrier  height  for  abstraction  because  the  meta¬ 
stable  intermediate  related  to  O2  [see  triangle  in  Fig.  4(a)] 
is  absent  [15].  To  investigate  whether  the  spin  selection 
rules  proposed  by  Ref.  [16]  are  necessary  to  account  for  the 
barrier,  we  additionally  calculated  the  singlet  EPES.  We 
find  minimal  differences  between  the  triplet  and  singlet 
EPESs,  which  indicates  a  very  small  energy  (—50  meV) 
required  to  change  the  spin  of  the  embedded  metal  cluster, 
and  hence,  of  the  entire  system.  Put  differently,  the  metal 
surface  easily  adopts,  within  reasonable  limits,  an  arbitrary 
spin  configuration,  thus,  substantially  decreasing  the  energy 
difference  between  different  (S2)  quantum  numbers  for  the 
entire  system.  Indeed,  at  room  temperature,  (S2)  of  an 
infinitely  extended  A1  surface  changes  under  the  influence 
of  thermal  fluctuations.  Lastly,  we  note  that  although  DFT 
fails  to  describe  dissociative  adsorption  of  02  on  Al(l  1 1), 
similar  calculations  for  seemingly  more  complex  transition 
metals  like  palladium  or  platinum  are  modeled  reasonably 
well  [33].  We  attribute  this  discrepancy  to  the  small  work 
function  of  Al(lll),  4.26  eV,  combined  with  the  delocal¬ 
ized  nature  of  its  valence  electrons,  which  greatly  facili¬ 
tates  charge  transfer  and,  coupled  with  DFT  XC  errors, 
leads  to  no  barrier  within  DFT.  The  higher  work  functions 
for  Pt(l  1 1)  of  5.93  eV  and  Pd(l  1 1)  of  5.6  eV  [34]  result  in 
an  energy  cost  for  charge  transfer  which  can  even  be  seen 
within  DFT. 

In  conclusion,  we  investigated  the  interaction  of  02  with 
Al(lll)  using  highly  accurate  embedded  correlated  wave 
function  methods  [23].  Unlike  conventional  KS  DFT 
calculations,  our  results  agree  with  experimental  observa¬ 
tions,  showing  finite  energy  barriers  at  all  high- symmetry 
points  of  the  surface,  for  both  parallel  and  perpendicular 
orientations  of  02.  The  barrier  arises  from  the  energetic 
cost  to  initiate  abrupt  charge  transfer  from  the  metal  surface 
to  the  molecule;  the  self-interaction  error  and  the  lack  of  the 
derivative  discontinuity  in  approximate  DFT  XC  cause 
spuriously  easy  charge  transfer  with  no  barrier  predicted. 
By  correctly  evaluating  charge  transfer  with  a  correlated 
wave  function  theory,  we  show  there  is  no  need  to  invoke 
nonadiabatic  surface  hopping  or  spin  selection  rules  to 
explain  the  barrier.  In  reality,  abrupt  spin  changes  at  the 
crossings  of  different  diabatic  02  spin  configuration  sur¬ 
faces  are  accommodated  by  small  spin  fluctuations  within 
the  metal  surface.  For  parallel  incidence,  the  potential 
energy  surface  points  to  direct  adsorption  at  two  adjacent 
fcc-hcp  hollow  sites  of  the  Al(l  1 1)  surface.  By  contrast,  we 
predict  abstraction  at  perpendicular  incidence:  our  results 
not  only  support  this  mechanism  as  dominant  pathway  at 
low  incident  energies  (as  proposed  by  STM  investigations), 
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but  also  provide  exquisite  detail  of  charge  and  spin  changes 
along  the  pathway.  Of  course,  conclusions  drawn  strictly 
from  PESs  may  miss  dynamical  details  such  as  steering 
events.  Future  quantum  dynamics  simulations  based  on 
our  EPES  may  confirm  the  proposed  mechanism. 
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